Checkerboard and stripe inhomogeneities in cuprates 
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We systematically investigate charge- ordering phases by means of a restricted and unrestricted 
Gutzwiller approximation to the single-band Hubbard model with nearest (t) and next-nearest neigh- 
bor hopping {t'). When \t' /t\ is small, as appropriate for La2-xSrxCu04, stripes are found, whereas 
in compounds with larger \t' /t\ (such as Ca2-xNaxCu02Cl2 and Bi2Sr2CaCu208+5) checkerboard 
structures are favored. In contrast to the linear doping dependence found for stripes the charge 
periodicity of checkerboard textures is locked to 4 unit cells over a wide doping range. In addition 
we find that checkerboard structures are favored at surfaces. 

PACS numbers: 71.28.-(-d,71.10.-w,71.45.1r,74.72.-h 
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The presence of charge and related spin inhomo- 
geneities in underdoped high-temperature supercon- 
ducting cuprates has received substantial experimen- 
tal support 0, S H Q. Less clear is the symmetry 
of the underlying textures, an issue which is strongly 
debatedHHlEla. 

Neutron scattering experiments on lanthanum 
cuprates (LCO) are usually interpreted in terms of 
one-dimensional modulations (stripes) where the two- 
dimensional symmetry of the scattering is explained with 
the presence of orthogonal stripe domains |5j. Indeed, 
a substantial anisotropy in the dynamic scattering 
between the two planar axis, as expected for fluctuating 
stripes, has been explicitely demonstrated in studies 
on untwinned samples of YBaaCugOy (YBCO)0, @. 
On the other hand, scanning tunneling microscopy 
(STM) in other classes of cuprates, iiamely the bilayer 
Bi2Sr2CaCu208+5 (Bi-2212)|i H (ni IjJ and the 
single-layer Ca2-xNaxCu02Cl2 (Na-CC0C)[l3 has 
revealed two-dimensional modulations (checkerboards). 
The peculiar characteristics of the Fermi surface of 
Na-CCOC ^3 is also in agreement with the expected 
features of a disordered checkerboard phase [l^. 

Charge-ordering (CO) phenomena thus seem to be 
common (and possibly generically present) in cuprates 
but the different symmetries found are rather puzzling 
and makes one wonder if one should not reinterpret the 
neutron scattering experiments in LCO and YBCO, in 
terms of two-dimensional textures also0. 

CO where predicted in cuprates before any experimen- 
tal detection [l5|. However, its cooperation or competi- 
tion with superconductivity is still debated. In this re- 
gard a possible direct relation between Tc and the in- 
tensity of incommensurate low-energy scattering (Id] de- 
mands for a deeper understanding of the physical mecha- 
nisms inducing CO and of the symmetry of the textures. 

In this work we address this issue on the basis of 
the one-band Hubbard Hamiltonian. It has been ar- 



gued that within this model, the ratio between nearest- 
(t) and next-nearest- {t') neighbor hopping t'/t < is 
the main electronic parameter characterizing the differ- 
ent cuprate families [tR Il8j|. We find a transition from 
one-dimensional to two-dimensional textures as \t' /t\ is 
increased. In addition, we find that the tendency to form 
checkerboard structures is increased at surfaces. Our re- 
sults provide a clue to interpret the contrasting experi- 
mental results regarding the symmetry of CO structures 
in different cuprate compounds. 

We solve the Hubbard model within a unrestricted 
Gutzwiller approximation (GA). The GA energy is min- 
imized with respect to unconstrained charge, spin and 
double occupancy distributions (for details see Ref. 19]). 
Generally we restrict to solutions without spin canting 
but we have checked the results lifting this restriction in 
several cases. In the case of stripes, once the unrestricted 
solution has been found, the symmetry of the solution has 
been imposed on a supcrcell and the problem has been 
solved in momentum space, allowing to obtain smooth 
curves of the energy as a function of doping. In the other 
cases, the curves are not smooth due to restricted sam- 
pling but, we estimate, that the finite size corrections to 
the energy of single points are negligible due to the large 
size of the cells (up to 16 x 16). 

The use of the present mean-field-like scheme to ad- 
dress the phenomenon of CO in cuprates requires fur- 
ther justification: i) Previous studies have shown that 
the GA applied to a three-band Hubbard model and in 
the present model captures the phenomenology of stripes 
in LCO cuprates [inj. In particular it accounts for the 
behavior of incommensurability and of the chemical po- 
tential as a function of doping. This requires an accu- 
rate evaluation of the small energy difference between 
different stripe solutions, ii) Response functions, due 
to fluctuations on top of these mean- field states, are in 
good agreement with spectroscopies!^ , both in the 
charge |2.1j and in the spin sector|2^. Hi) The accuracy 
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of the charge distributions can be directly tested com- 
paring with exact methods. We have found that the 
charge profile in the GA is indistinguishable from den- 
sity matrix renormalization group computations |25j re- 
vealing stripes on a 7 x 6 Hubbard cluster doped with 4 
holes, iv) Finally, the charge profile predicted within the 
GA in the three-band Hubbard modelji^l, has recently 
been found to be in excellent agreement with a charge 
sensitive measurement 0. 

All this establishes the GA as a very reliable technique 
to determine charge inhomogeneities in a strongly corre- 
lated system. Although, as any mean-field-like method, 
it cannot be trusted to determine the presence of long- 
range order (especially in the delicate spin-sector) it is 
expected to give a reasonably accurate estimate of the 
relative ground state energies of different textures, which 
is what we need here, and of the short-range correlations 
which determine it. 

For LCO we fix the Hubbard on-site repulsion U/t — 8 
and t'/t = —0.2. This parameter set and a very simi- 
lar one reproduce the spectrum of magnetic excitations 
both in the undoped and doped phases |2J|. The Hubbard 
interaction U and the nearest-neighbor hopping are not 
expected to depend significantly on the material given 
the similarity of superexchange interactions |0. Instead 
t' has been found to be quite sensitive to the cuprate 
family [13, E3- Photoemission studies in the insulating 
parent compounds have shown that there is substantial 
splitting among states with momentum (7r/2,7r/2) and 
(tt, 0) (we set the lattice constant a = 1). Detailed stud- 
ies of a single hole in one band models suggest that this 
splitting roughly scales with t' 18]. Taking t'/t — —0.2 
for LCOHI we estimate t' in Bi-2212 and CCOC by 
rescaling t' with the observed splitting In this way 

we obtain t' OAt for Bi-2212 and t' 0.5t for 

CCOC. The increase of \t' /t\ from LCO to Bi-2212 is in 
agreement with LDA results (l^l whereas for CCOC the 
present semiempirical estimate is significantly larger. 

Fig. n shows the most relevant textures found in this 
study, (a) Is the bond-centered (BC) stripe solution 
already reported in previous studies |2^. (b)-(d) are 
checkerboard textures. At low doping (n/j) and for large 
values of \t' /t\ we find the configuration (b) to be partic- 
ularly stable. It can be seen as a lattice of ferromagnetic 
(FM) polarons or "ferrous" and is refered to as small fer- 
ron checkerboard (SFC) . Each ferron consists of a 2 x 2 
plaquette with uniform magnetization and resides at the 
intersections of a crossed array of antiphase domain walls 
of the antiferromagnetic (AFM) order parameter. As 
doping increases the domain walls populate more uni- 
formly and evolve to a configuration that resembles a 
crossed array of stripes (c). At higher doping the size 
of the ferromagnetic islands increases, producing a large 
ferron checkerboard (LFC) structure (d). In addition, 
we have considered arrays of spin-polarons (not shown) 
consisting of single spin flips which locally bind the holes. 
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FIG. 1: (Color online) CO textures for U/t = 8. The length of 
arrows is proportional to the spin density whereas the squares 
are proportional to the hole density minus the hole density in 
the insulator, (a) Bond-Centered stripe texture with charge 
spacing d = 4, doping nu ~ 1/8, and t'/t = —0.2; (b) SF 
checkerboard with d = 4 at a; = 0.06, and t'/t — —0.5; (c) 
same as (b) with Uh = 1/8; (d) LP checkerboard with d — 4 
and t'/t = -0.5 at = 1/8 



Fourier transforming the charge and spin distribu- 
tions, we find that the most intense Bragg peaks for 
the = 4 checkerboards (SFC, LFC) are at = 
(±7r/2, 0), (0, ±7r/2) (charge) and at O, = (tt ± 7r/4, tt ± 
7r/4) (spin). As argued by Tranquada|3 the rotation be- 
tween Qc and Qs is a signature for a checkerboard phase 
since in case of stripes one would have Qc = (0±27r/(i, 0) 
and Qs = (tt ± 7r/d,0) or Qc = (0,0 ± 27r/d) and 
Qs = (0, TT ± Tr/d) depending on orientation. 

In Fig. 121 we report the values of the energy per site 
E for selected checkerboard (dashed lines) and stripe 
(solid lines) textures in case of t'/t = —0.2 (a) and 
t'/t — —0.5 (b). The main general conclusion one can 
draw from Fig. [5] is that a sizable value of the next- 
nearest-neighbor hopping t' stabilizes checkerboard so- 
lutions with respect to stripes, which are more stable at 
small \t' /t\. This is the most generic and relevant finding 
of this work and provides a rational on why stripes are 
observed in LCO materials (where t'/t —0.2), while 
STM experiments detect checkerboard structures in Na- 
CCOC and Bi-2212, where t'/t is known to be substan- 
tially larger 'T^. 

Fig. |2Ia) details on this issue by comparing the en- 
ergy of checkerboard with stripe solutions as a function 
of t'/t. Although most of the CO structures gain energy 
by increasing \t'\, it is clearly apparent that the checker- 
board solutions take a much greater advantage from this 
increase and become more stable at \t' /t\ > 0.35 ^ 0.45 
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depending on doping. 

Bond-centered checkerboards take advantage from the 
"ferronic" nodes of the charge texture, while site-centered 
checkerboards (not shown) lack this important feature 
and have much higher energies. This is a major difference 
with respect to stripe structures, where bond-centered 
and site-centered textures are nearly degenerate. An- 
other important difference lies in the doping dependence 
of the periodicity. The most stable stripe solution at low 
doping (cf. Fig. EJ has d ~ l/(2nft,) (with d « 4 - 10) 
providing an explanat ion |2 Oj for the well-known linear 
behavior between doping and incommensurability In 
contrast, the bond-centered checkerboards lose (mag- 
netic) energy in the domain walls and it is more conve- 
nient for them to adjust the charge in the FM plaquettes 
keeping the domain walls as short as possible. Indeed 
the checkerboard structures with d > 4 (not shown) are 
at higher energy in the whole doping range. This pro- 
vides an explanation for the fact that in Na-CCOC the 
charge periodicity is independent of dop ing and locked at 
d = 4: U\. For Bi-2212 d « 4 - 4.7 ES [13 has been 
reported in STM experiments. However, here the inter- 
pretation is more difficult due to the presence of lattice 
modulations, mesoscopic inhomogeneities, and less neat 
CO peaks. 

The SFC is particularly stable for large \t'/t\ at Uh = 
1/16 ~ 0.06 (cf. Fig. [Sis) when each ferron accommo- 
dates one added hole. This finding can be substantiated 
from an estimate of the electronic energy of an isolated 
2x2 plaquette. One obtains the eigenvalue structure 
—2t — t';t']t';2t — t' with the last level unoccupied due 
to the presence of the hole. This yields the energy per 
site E — Eafm — Ja — t{l — t'/2t)/8 where —Ja is the 
energy per bond of the AFM phase. Additional carriers 
have to overcome a large gap {2t — 2t') leading to a rapid 
rise of the energy (cf. Fig.|2Il. At rih = 1/16 other con- 
figurations take a smaller advantage from t' < while 
the magnetic energy cost is higher in the checkerboard. 
Therefore a substantial value of \t' /t\ is needed to stabi- 
lize checkerboards with respect to stripes and polarons. 

At higher doping and for large \t'\ it pays to break 
more AFM bonds to have larger ferromagnetic islands 
to accommodate the holes. In this sense the SFC/LFC 
states locally reflect the tendency of the extended Hub- 
bard model towards ferromagnetism for large \t'\ |2^E^ 
and can be seen as inhomogeneous precursors of the uni- 
form ferromagnetic phase, although they do not have a 
net ferromagnetic moment. 

At intermediate dopings Uh ~ 0.1 and t'/t — —0.5 
a phase separated solution between the SFC and LFC 
has lower energy than c? = 4 stripes. Ignoring the long- 
range Coulomb interaction the energy of this solution is 
given by the Maxwell construction [dashed black line in 
Fig. Elb)] . This solution will be frustrated in the pres- 
ence of the long-range Coulomb interaction giving rise to 
a mesoscopic phase separation 28J . However, since the 
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FIG. 2: (Color online) Energy per site as a function of doping 
for the different textures studied and for t'/t = —0.2 (a) and 
t'/t = —0.5 (b). For clarity we subtracted the energy of the 
AFM solution and the line fionh with fio = —1.6t. Different 
choices of no correspond to different choices of the origin of 
the energy of the single particle states and do not change the 
physics. In each curve for the stripes the charge periodicity 
perpendicular to the stripe is fixed and takes the values (from 
left to right) d = 10, 8, 6, 5, 4, 3. In panel (b) we also show 
the energies allowing for spin canting and AFM-SFC mixing 
(filled circles) and Maxwell construction (MC) for AFM-SFC 
phase separation and for SFC and LFC phase separation. 

solutions have the same symmetry they can be mixed at 
a quite small length scale which implies a low Coulomb 
cost. 

At doping rifi < 1/16 the SFC is expected to phase 
separate with the AFM solution. MC for this case is also 
shown in Fig. [2Ib) while the corresponding charge- and 
spin distribution for rih — 0.047 is shown in Fig. |2|^b). 
In this case the two phases are quite different which 
can imply a high energy cost for mesoscopic phase 
separation|28j since a large surface energy may be ex- 
pected at the interface. We have found, on the contrary, 
that allowing for spin canting one finds a solution with 
negligible surface energy cost as shown by the filled cir- 
cle at Hh = 0.047 in Fig. H^b). Indeed, although the 
solutions have a substantial amount of interface due to 
the finiteness of the cluster [Fig. El^b)] , the energy is very 
close to the MC line, for which, the surface energy cost is 
assumed to be zero. This is because the surface energy is 
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FIG. 3: (Color online) (a) Energy per site vs. ~t' /t for U /t — 
8,10 for stripes and checkerboard at Uh = 0.0625,0.187. 
(b)Mixed AFM-SFC solution at nu = 0.047 and allowing for 
spin canting. 



dominated by magnetic effects, but spin twisting allows 
for a smooth "flipping" of the magnetic order parameter 
from the SFC to the pure AFM in accord with standard 
arguments for an order parameter which breaks a con- 
tinuous symmetry 01 . We also see from Fig. I^b) that 
for the pure checkerboard solution at Uh = 1/16 the spin 
canting energy gain is negligible. 

Mesoscopic textures like the one shown in Fig. Efb) 
will also be influenced by disorder effects which will ran- 
domly pin clusters of one or the other phase according to 
their charge. It is clear from Fig.Ofb) that this will also 
induce substantial disorder on the spin degrees of free- 
dom providing a natural mechanism for spin-glass effects 
often seen in underdoped cuprates. 

As in previous studies we find that with a ratio of \t' /t\ 
appropriate for LCO the more stable inhomogeneities are 
stripes. As \t' /t\ is increased checkerboard structures are 
favored. Bi-2212 and Na-CCOC are estimated to lie close 
to the transition 30] providing an explanation for the puz- 
zling difference in symmetry observed. The limited ac- 
curacy of specific parameter estimates and of the critical 
\t'/t\ for checkerboard structures does not allow to ex- 
clude stripes in the bulk of these materials. Quantum 
fluctuations may also blur the difference among different 
configurations close to the transition. 

So far checkerboard textures have only been detected 
by surface sensitive probes but are not yet established 
as a bulk phenomenon. Since the Hubbard U interac- 
tion is screened from its atomic to the solid-state value 
by the polarization of the environment '3 1] a substantial 
increase of U is expected to occur at surfaces. Increasing 
the value to U/t = 10 we find [cf. Fig|3ta)] that checker- 
board structures are more favored with respect to stripes 
and the critical value of \t' /t\ is decreased by about 15%. 
It is thus possible that checkerboard structures are an ex- 
ample of electronic surface reconstruction and occur only 
at surfaces 32]. 

In conclusion, we have shown that the next-nearest 
neighbor hopping plays a dominant role in determin- 
ing the symmetry of CO textures in cuprates. Whereas 



for a ratio t' /t = —0.2 as appropriate for LCO we find 
stripes as the most stable inhomogeneities, a crossover to 
checkerboard solutions occurs for large but still realistic 
values of \t'/t\. Checkerboards behave quite differently 
as a function of doping with respect to stripes in that 
the periodicity of the modulation is locked as observed 
experimentally^^. The competition between FM clus- 
ters and antiphase domains rules the relative stability of 
the various textures. Spin canting produces a negligible 
energy gain for ordered structures but it is fundamental 
for nanoscopically mixed phases. 

A change of symmetry from stripes in the bulk of LCO 
and YBCO to checkerboards at the surface of Bi-2212 
and Na-CCOC appears quite naturally, thereby explain- 
ing the puzzling lack of universality found for the sym- 
metry of CO. On the other hand CO itself appears to be 
quite ubiquitous in cuprates. 
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